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In  Situ  Optimal  Reshading  of  Arrays  with  Failed  Elements 


ALGORITHM  DOCUMENTATION  PACKAGE  j 

This  document  assembles  under  one  cover  information  on  a  NUSC-developed  algorithm  which 

computes  optimal  shading  weights  for  discrete  elements  (sensors)  in  linear  acoustic  arrays.  The 

algorithm  has  been  found  especially  useful  when  elements  fail  and  array  reshading  is  required  in  ! 

situ.  The  main  attractions  of  the  algorithm  are  that  it  loads  easily  on  Hewlett-Packard  microcom-  *j 

puters,  and  that  it  runs  fast  enough  and  is  accurate  enough  to  suit  most  sea  trial  and  engineering 

development  applications.  Continuing  requests  for  this  information  since  an  invited  paper  first 

appeared  in  the  IEEE  Journal  of  Oceanic  Engineering  in  January  1987  motivated  the  publication  of 

this  documentation  package. 

• 

The  information  included  here  is  in  hard  copy  and  floppy  disk  form:  the  IEEE  paper.  In  Situ  Optimal 
Reshading  of  Arrays  with  Failed  Elements,  is  reprinted  in  Appendix  A;  a  program  listing  of  the 
application-specific  driver  routine  is  given  in  Appendix  B;  and  a  3  Vz  inch  floppy  disk,  containing  the 
program  “Reshade”  which  runs  on  any  Hewlett-Packard  Series  200  or  300  micrt  -omputei,  is  pocketed 
in  Appendix  C.  , 


GENERAL  APPIJ  CATION 

When  a  linear  array  of  discrete  acoustic  elements  is  subjected  to  the  rigors  of  the  ocean  environment, 
individual  elements  can  fail.  Element  failures  are  usually  characterized  by  noisy  channels,  or 
intermittent  responses,  or  no  response  at  all.  Depending  on  the  number  of  failed  elements  and  their 
specific  locations  within  the  array,  sidelobe  levels  of  the  array  wavenumber  (it)  response  can  rise 
significantly  to  degrade  array  performance.  Because  element  weighting  values  determine  array 
wavenumber  response,  weights  that  are  optimal  for  a  fully  populated  array  have  to  be  recalculated 
when  elements  fail.  The  optimal  reshading  (reweighting)  algorithm  described  here  can  be  applied 
in  situ  to  compute  weighting  values  that  can  reduce  sidelobe  levels  to  approximately  the  original 
design  specification.  In  fact,  in  the  more  common  situations  where  “a  few”  elements  fail,  optimal 
reshading  does  regain  original  sidelobe  levels.  Where  large  numbers  of  elements  fail,  optimal 
reshading  is  still  possible  but  may  be  of  limited  use. 

The  original  approach  for  optimal  reshading  of  a  linear  array  was  proposed  by  Streit  and  Nuttall  in 
1982  (see  reference  1,  Appendix  A).  At  that  time,  the  algorithm  was  run  on  a  VAX  11/780  and 
required  hours  of  computation  time  and  large  amounts  of  mass  storage  for  rudimentary  element 
failure  problems. 

The  1987  reshading  algorithm  incorporates  several  algorithmic  improvements  that  exploit  the  special 
structure  of  the  underlying  linear  programming  problem  to  reduce  time  and  storage  requirements 
by  orders  of  magnitude.  The  current  algorithm  is  still  based  on  the  original  theory,  but  is  now  fast 
enough  and  small  enough  to  execute  successfully  in  minutes  instead  of  hours  in  the  application 
environment.  Execution  time  for  a  50-element  array  is  typically  about  10  minutes.  Derivation  of  the 
optimal  reshading  algorithm  and  its  implementation  are  given  in  the  references  of  the  paper  reprinted 
in  Appendix  A;  examples  of  array  reshading  are  given  in  the  paper  itself. 
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TO  RUN  “RESHADE”-. 

•  Insert  the  program  disk  from  Appendix  B  into  device/drive. 

•  Type  the  command  string  LOAD  “ RESHADE : (device  specifier)* 

•  Press  Enter 

•  When  the  program  is  loaded,  press  Run. 

•  Follow  the  prompts. 


PROGRAM  NOTES 

“Reshade”  comprises  a  driver  routine  in  HP  BASIC  which  sets  up  the  necessary  variables  to  be 
optimized  and  a  generic  optimization  routine.  The  driver  is  listed  in  lines  1  through  481  of  the 
program — the  printout  is  contained  in  Appendix  B. 

The  driver  included  in  “Reshade”  applies  to  a  linear  array  of  acoustic  elements,  some  of  which  may 
have  failed  during  the  course  of  a  sea  trial  or  similar  event.  Even  with  the  array  intact,  “Reshade” 
allows  the  user  to  minimize  the  sidelobe  levels  of  the  array  beamformer  output,  given  a  certain 
mainlobe  width.  If  the  minimum  sidelobe  levels  remain  n>o  high,  it  is  possible  to  alter  the  mainlobe 
beamwidth  to  reduce  sidelobe  levels.  Note  that  the  weights  on  each  element  can  be  set  up  as 
non-negative,  if  desired. 

The  program  prompts  require  user  inputs,  not  all  of  which  are  self-explanatory.  For  each  user  input, 
values  in  (parentheses)  are  those  allowed,  and  values  in  [brackets]  are  the  defaults.  The  maximum 
allowable  total  number  of  array  elements  is  50;  the  minimum  is  three.  The  computation  time  for  a 
50-element  array  is  approximately  10  minutes,  while  a  10-element  array  nuts  in  less  than  one  minute. 

The  algorithm  is  applicable  to  both  equispaced  or  aperiodically  spaced  linear  arrays.  In  the 
equispaced  arrangement,  the  wavenumbers  kO  and  kl — which  delimit  the  region  in  which  the 
minimization  is  performed — are  calculated  automatically  from  the  desired  sidelobe  level.  The  final 
sidelobe  level  depends  upon  the  number  of  failed  elements  in  the  array  and  their  location.  In  this 
case,  only  the  inter-element  spacing  must  be  specified. 

In  the  aperiodically  spaced  arrangement,  every  element’s  position  referenced  to  the  forward  end  of 
the  array  must  be  specified.  If  elements  have  failed  (or  are  missing),  they  are  treated  as  if  they  do 
not  exist.  The  wavenumbers  kO  and  kl  are  not  calculated  automatically  for  an  aperiodic  array,  and 
must  be  entered  manually  in  units  of  radians/metet 
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If  unsatisfactory  sidelobe  levels  are  still  present  after  running  “Reshade”  for  an  equispaced  arrange¬ 
ment,  kO  can  be  increased  to  provide  a  larger  beamwidth,  thus  reducing  sidelobe  levels.  For  an 
aperiodic  anay,  kO  and  kl  can  be  altered  manually  to  reduce  the  sidelobe  level  in  the  region  of 
interest. 

Resultant  weights  can  be  stored  in  a  data  file  in  the  following  format: 

•  Equispaced  element  arrangement — total  number  of  elements,  followed  by  the  inter-element 
spacing,  followed  by  array  weights. 

•  Aperiodic  element  arrangement — total  number  of  elements,  followed  by  each  element’s 
position,  followed  by  array  weights. 


PROGRAM  EXTENSIONS  AND  IMPROVEMENTS 

“Reshade”  and  its  associated  algorithm  have  established  the  validity  of  in  situ  computation  of  linear 
acoustic  array  optimal  shading  weights.  Virtually  no  sea  trial  is  conducted  today  without  reshading 
to  compensate  for  failed  elements.  Extensions  to  larger  linear  array  problems  are  potentially  useful. 
Improvements  and  modifications  to  the  original  source  code  are  possible  in  the  light  of  recent 
advances  in  signal  processing  hardware,  and  are  needed  to  obtain  reasonable  computation  times  for 
these  larger  arrays.  With  the  advent  of  single-board  array  processors,  the  beam  pattern  computations 
done  (implicitly)  in  each  iteration  in  the  generic  optimization  model  (KAPROX)  may  be  performed 
more  quickly  and  accurately  using  a  floating  point  FFT  This  is  but  one  example  of  software 
modifications  which  will  enhance  the  performance  o»  “Reshade.” 

The  generic  nature  of  the  optimization  routine  lends  itself  to  the  solution  of  more  general  array 
problems.  These  arrays  may  be  multiline,  planar;  or  three-dimensional  with  arbitrary  geometry.  Each 
geometry,  however,  will  require  a  specific  driver  routine  to  set  up  the  problem  to  be  optimized.  In 
general,  the  drivers  would  need  the  capability  to  address  complex  weights,  allocate  enough  memory 
for  computations,  and  to  take  into  account  any  application-specific  constraints  imposed  on  the 
optimization  problem.  Additional  constraints  can  be  useful;  for  instance,  constraints  can  sometimes 
be  used  in  active  arrays  to  control  adverse  effects  of  acoustic  coupling  between  the  anay  elements. 
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In  Situ  Optimal  Reshading  of  Arrays  with  Failed  " 

Elements 

MICHAEL  S.  SHERRILL  and  ROY  L  STREIT,  senior  member,  jefe 

{Invited  Paper )  i 


Abut  act— An  itgorilba  is  presented  which  computes  optimal  weights 
for  arbitrary  linear  arrays.  The  application  of  this  algorithm  to  in  situ 
optimal  reshading  of  arrays  with  failed  elements  is  discussed.  It  is  shown 
that  optimal  reshading  can  often  regain  the  original  sidelobe  level  by 
slightly  increasing  the  mainlobe  beamwidth.  Three  examples  are  pre¬ 
sented  to  illustrate  the  algorithm's  effectiveness.  Hardware  and  software 
issues  are  discussed.  Execution  time  for  a  25-elcmen(  array  is  typically 
between  1  and  2  min  on  an  HP9836C  microcomputer. 

I  Introduction 

linear  array  of  discrete  elements  (sensors)  often  experi¬ 
ences  element  failures  in  situ.  These  failures  can 
significantly  increase  the  sidelobe  levels  of  the  array 
wavenumber  response,  depending  on  how  many  elements  fail 
and  where  the  elements  are  located  within  the  array.  We 
discuss  here  an  optimal  reshading  (reweighting)  algorithm 
which  can  be  applied  in  situ  to  reduce  the  sidelobe  levels  to  the 
original  design  level.  In  many  common  element-failure  situa¬ 
tions,  optimal  reshading  can  regain  the  original  sidelobe  level 
by  slightly  increasing  the  mainlobe  beamwidth.  In  arrays 
which  experience  significant  element  failures,  optimal  reshad¬ 
ing  is  still  possible,  but  may  be  of  limited  use.  Three  examples 
given  below  demonstrate  a  few  of  the  possibilities. 

An  algorithm  for  optimal  reshading  was  first  proposed  in  ( I ) 
by  Streit  and  Nuttall.  Their  algorithm  utilized  the  general- 
purpose  subroutine  (2|  to  solve  a  specially  structured  "linear 
programming"  problem.  Unfortunately,  their  algorithm  re¬ 
quired  hours  of  computation  time  and  large  amounts  of 
computer  storage  on  a  minicomputer  (the  VAX  11/780)  to 
optimally  reshade  a  50-element  array  with  five  failed  ele¬ 
ments,  Consequently,  their  algorithm  is  not  useful  for  in  situ 
optimal  reshading. 

The  shading  algorithm  proposed  here  differs  from  Streit  and 
Nuttall's  primarily  in  that  we  solve  their  linear  programming 
piobiem  using  a  new  general-purpose  subioutine  (5J,  i , 
herein  referred  to  as  Algorithm  635.  Algorithm  635  uses  the 
special  structure  of  the  linear  programming  problem  to  reduce 
time  and  storage  requirements  by  orders  of  magnitude. 
Algorithm  635  can  be  incorporated  easily  in  Streit  and 
NutUll'i  original  approach.  A  significant  algorithmic  im¬ 
provement  was  discovered  in  the  course  of  this  study  and  is 
described  below.  The  resulting  shading  algorithm  is  fast 
enough  and  small  enough  to  execute  successfully  on  micro- 
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computers  (such  as  the  HP9836C  used  here)  in  only  a  few 
minutes.  Typical  execution  time  for  a  25-element  array  is 
under  2  min;  for  a  50-element  array,  execution  time  is 
typically  under  10  min  The  current  algorithm,  and  the 
HF9836C  with  ns  inherent  transportability,  comprise  an 
effective  system  for  optimal  reshading  in  situ 

II  Optimal  Array  Shading 

The  wavenumber  response  of  a  linear  array  composed  of  A’ 
discrete  omnidirectional  elements  located  at  arbitrary  fixed 
positions  x„  is  given  by 

V 

r<*)=  £  wr  exp \-ikxn\  (1) 

where  w„  are  the  element  weights  and  the  independent  variable 
k  denotes  wavenumber  in  radians  per  unit  length  The  element 
weights  are  required  to  be  real,  but  this  entails  no  loss  of 
generality  (see  below  in  Section  III)  Also,  from  (1).  T(-k ) 
=  T*(k)  for  real  weights  (asterisk  denotes  conjugation),  so  it 
is  unnecessary  to  consider  negative  values  for  k  and  wc 
confine  our  attention  to  nonnegative  k 

The  array  response  as  a  function  of  k  can  be  considered  to 
be  composed  of  a  mainlobe  beamwidth  and  a  sidelobe  region 
The  objective  of  the  optimization  process  is  to  make  j  T(k)j  as 
small  as  possible  on  the  user-specified  sidelobe  interval.  Array 
weights  which  achieve  this  objective  are  said  to  be  optimal 
The  optimization  process  usually  produces  cquivalued  side- 
lobes  in  the  sidelobe  region. 

Weights  that  arc  optimal  for  a  full  array  do  not  remain 
optimal  after  the  array  experiences  element  failures  To 
oartially  compensate  for  failed  elements,  the  array  is  optimally 
reshaded  by  undertaking  the  optimization  process  again  and 
incorporating  knowledge  of  which  elements  have  failed  As 
the  examples  Vlow  will  •n*  effectiveness  of  this 

strategy  depends  upon  how  many  elements  have  failed  and  the 
location  of  these  elements  in  the  array. 

The  sidelobe  interval  is  defined  differently  depending  on  the 
interelement  spacing  of  the  array  For  an  array  with  periodi¬ 
cally  spaced  elements  and  no  failures,  the  sidelobe  interval  is 
defined  to  be  [AZ0,  (2i/D)~  A'0|,  where  K0  is  calculated  from 
the  desired  sidelobe  level  ano  the  number  N  ..f  array 
elements  1  D  is  ihe  physical  distance  from  sensor  to  sensor 

1  For  *n  A-elemem  trray  and  -  r-dB  peak  sidelobes,  we  have  K0  =  (2/D) 
arccos  (l/Z0>  where  2 Z.  =  |r  -  (r1  -  l)1'2)''*'  +  (r  +  (r1  -  t)*'JJ''", 
r  -  to"®,  and  M  =  /V  -  1  The  interelement  spacing  D  is  assumed  to  be 
half  of  the  so-called  design  wavenumber,  and  /V  is  ihe  number  of  array 
elements  before  failures 
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Furthermore,  the  minimization  interval  can  be  reduced  f«>  |  A, . 
*■  D\,  since  the  response  of  this  array  is  symmetric  about  A 
x  D  A'ti  is  typically  the  point  on  the  mainiobe  response  which 
is  equal  m  magnitude  to  the  peak  of  the  sidelobes.  hut  this  is 
not  always  true  for  seriously  degraded  and  or  aperiodic  arravs 
(see  Example  3  below)  For  arrays  with  apenodically  spaced 
elements,  the  sidelobe  interval.  denoted  by  j  A ...  A  }.  must  be 
chosen  by  inspection  of  a  nonoptimai  beam  pattern  or  some 
other  means  I  T{k)  must  be  minimized  over  the  full  1  A  ,  A  | 
range  since,  in  general,  an  aperiodic  array  response  is  not 
symmetric  about  any  wavenumber  other  than  A  =  0  The 
ability  to  specify  arbitrary  AVi  and  K<:  is  particular!)  useful  for 
those  applications  involving  apenodically  spaced  elements 
because  lower  sidelobe  levels  may  be  obtained  by  looking  at 
different  minimization  regions 

The  optimization  process  deals  with  element  failures  in  an 
array  in  the  tollowing  way 

Step  1  Maintain  mamlobe  beam  width  and  permit  the 
sidelobe  levels  to  rise 

Step  2  Regain,  if  possible,  the  original  sidelobe  level  b\ 
broadening  the  mamlobe 

Broadening  the  mainiobe  by  increasing  A  i stop  2i  is  per¬ 
formed  only  if  the  sidelobe  level,  even  after  optima!  reshad 
mg.  has  risen  to  an  unacceptable  value  because  of  element 
failures  Thus  step  I  is  normal  algorithmic  procedure,  and  step 
2  requires  some  iteration  in  specifying  A.  and  or  A  because  a 
compromise  has  to  be  made  between  the  mainiobe  beamwidth 
and  the  level  of  the  sidelobes 

The  solution  of  the  array  problem  in  the  original  formula 
tion  1 1 1  is  mathematically  equivalent  to  solving  an  overdeter 
mined  system  of  complex  linear  equations  Unacceptably  high 
sidelobes  result  if  this  system  is  solved  in  the  usual  leas! 
squares  sense,  so  it  is  necessary  to  solve  the  system  so  that  the 
magnitude  of  the  maximum  residual  error  is  minimized  There 
now  exists  ( 3 ]  an  efficient  algorithm  and  corresponding 
FORTRAN  axle  [4|  tor  solving  problems  of  this  son  to  high 
accuracy 

To  obtain  the  beamformer  equation  in  an  appropriate  format 
to  utilize  this  algorithm,  we  normalize  the  peak  response  of 
T{k)  so  that  7(0)  --  1.  This  gives 

2  w*  =  1  (2) 
*  »  I 

We  solve  (2)  for  the  Vth  weight  and  substitute  in  (1)  to 
obtain 

iV  -  I 

T(k)  =  cxp  |  -  rfar.v)  +  2  *5,  [exp  f  "  tfar*)  -  exp  (  -  /Aav)1 

(3) 

By  sampling  T{k)  at  the  M  equispaced  points 
„  [*i  -  A0J 

km-  A0+ - (m  -  1),  m  ~  1 ,  •  •  ,  M  (4) 

M  -  I 

we  can  write  the  problem  of  minimizing  the  peak  sidelobe 


level  of  the  array  response  as 

nun  max  V  u,„..  n  „  i  ^ } 

where  the  complex  numbers  Jm  and  a...r  are  defined  K 

{..  exp  |  ik„  i  I 

d„.r  exp  S  /A.,.  1  exp  j  ;4.„.  j 

The  problem  (5)  is  precisely  the  form  necessary  for  application 
of  Algorithm  635  For  theoretical  details  of  this  algorithm,  the 
interested  reader  is  referred  to  |3j 

Sometimes  a  few  of  (he  optimum  weights  for  array-,  with 
failed  elements  are  observed  to  be  negative,  pjnivularjv  those 
on  the  end  elements  I?  the  weights  are  applied  in  hardware, 
providing  a  180°  phase  factor  on  the  element  output  mav  not 
be  desirable  or  possible  However.  Algorithm  allows  the 
selection  of  all  nonnegative  weights,  th.:-  is  implemented  by 
the  addition  of  constraints  to  1 5 »  Usually  hut  not  always,  art 
element  is  ze rived  if  if  would  have  had  a  negative  weigh!  From 
(2)  it  follows  that,  if  all  the  element  weight  values  are  required 
to  be  positive,  they  must  be  between  0  and  I  The  requirement 
that  weights  w  .  .  w be  between  0  and  I  can  be  written 

mathematically  as 

l  1 

n  1  .  ,  \  I  { "  i 

Algorithm  635  requires  these  .V  l  constraints  Algorithm 
635  can  also  incorporate  any  number  of  general  constraints  of 
the  form 

V  6„.„v i„  <•„.  <  d„  .  rn  1.2.  .  /  i  S  > 

where  c„,  and  d,„  are  constants  The  require.. tent  that  * also 
be  nonnegdtive  gives 


or 

I  I 

s  »•„--<•  (9) 

4-*  S  S 

which  is  clearly  a  special  case  of  the  general  constraints  (8i 
III  Algorithm  Improvimlnis 

Several  changes  to  the  algorithm  presented  in  [  1 1  enable 
significant  reduction  in  the  need  for  computational  intensity 
Lewis  and  Strctt  [5]  !.«>c  proved  fh3»  for  a  general  line  arrav 
shaded  so  that  it  has  optimal  sidelobe  levels  when  steered 
through  the  same  number  of  degrees  either  side  of  broadside, 
there  exists  a  set  of  optimal  weights  that  are  real  Thus 
complex  weights  do  not  need  to  be  considered  This  fact 
allows  an  approximate  eight  fold  reduction  in  computation 
time  and  a  two-fold  reduction  in  storage  requirements 


\ 
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It  is  clour  that  the  50-element  example  run  in  Streit  and 
N. Lilt Ji  1  |IJ  was  significantly  oversampled  in  wavenumber 
Their  beam  pattern  can  he*  reproduced  « nh  a  tour- fold 
reduction  in  the  sampling  of  T\k) ' see  Kxample  2  below  i.  and 
this  m  no  was  detracts  from  the  practical  application  of  the 
algorithm  A  significant  reduction  in  compulation  time  is 
realized  b>  decreasing  the  number  \f  of  beam  pattern  samples 
in  i4) 

A  significant  algorithmic  modification  made  to  \inthm 
635  further  decreases  computation  time  We  h  v*'  unrled  this 
modification  "fast  costing”  and  it  is  an  ,*i.t«»rtant  step  in 
making  the  algorithm  feasible  on  microcomputers  such  as  the 
HP^lObC  In  order  to  describe  »'ii»  modification  properly, 
k'tne  familiarity  with  the  smif.ex  method  of  linear  program 
mine  and  reference  |3]  is  ssumed 

Algorithm  635  car  Voken  into  two  fundamental  compu 
rational  operations  called  “costing''  and  “pivoting  “Cost 
mg"  determines  the  so  called  minimum  reduced  cost  coeffi 
^ lent  and  requires  2S\t  multiplications,  where  .V  i>  the 
number  of  discrete  array  elements  and  M  is  the  number  of 
sampies  taken  of  the  beam  pattern  Pivoting’  is  a  basis 
update  and  requires  \  *  real  muinplication.s  It  is  Jear  that  the 
speed  ot  the  algorithm  is  mtinutdv  related  to  the  number  M  of 
sample'  taken  of  the  beam  pattern,  as  well  as  the  number  ,\  <»} 
discrete  array  elements  Since  M  is  larger  than  .V.  costing" 
require'  more  multiplications  ihun  “pivoting 

Costing’  in  the  linear  arras  application  means  that,  in  each 
'implex  iteration,  the  discreti/cd  absolute  value  ’  of  every 
sidelobe  sample  of  the  wavenumber  response  function  Tt  km). 
tv  I,  .  M,  »s  computed  to  deteiinine  the  “minimum 
reduced  sost  coefficient”  ot  the  current  “basic  feasible 
solution  ’  B>  proceeding  through  a  finite  sequence  ot  such 
'basic  feasible  solutions,  we  arrive  at  the  solution  of  the 
“discretized  problem  “  As  shown  in  {3j.  this  implies  that  the 
computed  optimal  wavenumber  response  function  can  have 
sidelobe  levels  that  are  theoretically  at  most  0  04  dB  higher 
than  the  true  optimum  sidelobe  level  :  Fast  costing"  refers 
simply  to  the  fact  that  we  first  determine  which  of  the  sidelobe 
samples  T{km),  m  =  1.  •  M.  has  the  largest  true  absolute 
value,  and  then  compute  the  “discretized  absolute  value"  of 
this  one  complex  number  Therefore,  only  one  “discretized 
absolute  value"  calculation  is  performed  in  each  simplex 
iteration  instead  of  M  such  calculations  The  resulting 
reduction  in  computational  effort  is  significant  in  microcom¬ 
puting  environments  The  drawback  is  thal  the  use  of  "fast 
costing"  prevents  the  simplex  algorithm  from  converging  to  a 
solution  of  the  "discretized  problem  "  Fortunately,  however, 
it  can  be  proved  that  we  must  approximate  the  solution  in  a 
well  defined  sense  In  the  linear  array  application,  "fast 
costing"  results  in  the  computed  optimum  beam  pattern 
having  sidelobe  levels  that  arc  theoretically  at  most  0  OH  dB 
higher  than  the  true  optimum  level  •  This  is  a  small  price  to 
pay  for  major  execution  time  improvements 

:  The  theoretical  error  of  at  most  0  04  dB  iv  derived  by  uking  20  log.0  (tec 
•  r  pi),  where  p  ■-  32  The  term  sec  f  9/p)  n  (he  error  bound  discussed  tn 

IB 

fast  costing  squares  the  error  bound,  giving  sec;  <  *'p).  or  0  Ok  dB  when 
P  32 


IV  AlOOkMHM  Istut-su  M  A'lTiN  |MK  f\  Stn  f 
An  algorithm  must  be  reliable.  e.j'\  !o  use  and  fjst  when 
executing  on  portable  microcomputers.  to  he  useful  tor  in  vitu 
application  The  following  section  details  the  most  important 
hardware  and  .software  issues  addressed  to  enable  in  \itu 
optimal  reshading  of  arrays  with  faded  elements 

The  algorithm  has  been  coded  ir  BASIC  and  is  comprised  <4 
Algorithm  635  and  an  array  processing  driver  program 
Algorithm  MS  solves  the  linear  program  tor  a  set  ot  optimal 
weights,  given  data  supplied  by  the  driver  program  I  he 
driver  performs  the  initial  setup  based  on  several  user  enpu/s 
and  provides  ail  piogum  output 

The  driver  program  may  he  used  with  linear  arrays  having 
either  periodic  or  jperiodically  spaced  elements  Program 
output  consists  of  a  graph  of  the  optima!  beam  pattern,  a  graph 
of  the  optimal  normalized  element  weights,  and  several 
parameters  pertinent  to  the  specific  problem  Provision  is 
made  for  storing  the  weights  in  a  separate  data  file  tor  p^s'sHc 
use  with  digital  heamformers 

A  Hewlett  Packard  »HP»  specific  software  m«*difuation  was 
made  by  setting  up  the  input  data  arravs  (equation  Coi  in 
buffers  so  that  they  are  accessible  for  a  one  dimensional 
multiplv  For  large  array  dimensions,  indexing  a  doubly 
subscripted  data  array  and  performing  a  dot  product  takes 
more  time  on  the  Hfx)H  *6C  than  reading  in  a  data  arras  from  j 
buffer,  doing  a  MAI  multiply,  and  performing  a  summation 
(A  MAT  multiply  is  simply  an  element  hv  element  multiplv  ot 
two  equally  dimensioned  data  arrays  i  However,  this  pone 
dure  is  more  time  consuming  when  the  input  data  arravs  are 
very  small  tie.  the  number  of  elements  in  the  line  array  i- 
small)  The  break  even  point  »vcurs  at  around  12  or  I* 
elements,  so  it  was  decided  to  insorporate  this  speed  enhance 
menl  for  the  longer  running  larger  line  arrays  and  trade  oil 
some  speed  reduction  on  the  smaller  line  arrays 

To  obtain  fast  execution  times  for  in  situ  applications,  we 
use  one  hardware  speed  enhancement,  a  12  5 -MHz  fast  CPI 
card  with  16  kbytes  of  cache  memory  This  hardwjre 
supplement  is  available  from  HP  for  use  on  the  HP9H36C 
Cache  memory  is  fast  memory  resident  on  the  CPC  card  for 
quick  instruction  acquisition  The  use  of  the  fast  CPI  board 
rather  than  the  8  MHz  clock  present  in  the  standard  computer 
configuration  results  m  an  approximate  factor -of- two  increase 
in  observed  speed 

The  complete  program  is  precompiled  by  use  of  software 
and  a  floating  point  math  card  available  from  the  INFOTFK 
company  Prccompilation  reduces  mo*  computational  por¬ 
tions  of  the  BASIC  code  to  machine  language,  giving  an 
additional  three  fold  reduction  in  computation  time  It  is  also 
desirable  to  upgrade  the  operating  system  for  the  HP  to  its 
latest  revision  All  work  on  these  problems  was  run  using  the 
BASIC  3.0  operating  system  and  the  hardware  supplements 
noted  above 

Computation  time  is  defined  as  time  spent  in  Algorithm  635 
and  does  not  include  the  small  amount  of  set  up  time  required 
by  the  driver  program  Computation  times  are  for  the 
compiled  BASIC  program  run  on  the  HP9836C  with  the 
special  hardware  additions  mentioned  above 

The  program  described  here  needs  just  over  303  kbytes  of 
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internal  memory  in  addition  to  the  memory  required  by  the 
operating  system  to  execute  on  the  HP9836C  Phis  is  the 
amount  of  space  required  by  fixing  the  maximum  array  size  at 
•V  -  50,  and  allowing  at  most  \f  -  256  beam  pattern 

samples  Users  can  change  dimensions  to  suit  their  specific 
needs,  but  storage  requirements  presently  are  directly  propor 
tional  to  the  product  A :\1  Even  tor  a  much  larger  number  or 
line  array  elements,  it  is  unlikely  that  memory  restrictions 
would  prove  to  be  a  problem  on  the  HP9836C  since  extra 
memory  boards  of  l  Mbyte  each  are  readily  available. 

Ongoing  modifications  should  further  enhance  the  capabil 
ity  and  speed  of  the  BASIC  algorithm  and  driver  The  addition 
of  the  ability  to  handle  directional  sensors  is  both  useful  and 
straightforward  to  implement  Execution  of  identical  code  on 
the  new  HP  300  series  computers,  which  have  a  16  6  MHz 
clock  rate,  should  further  reduce  the  computation  time 
Compulation  times  on  the  order  of  5  min  for  a  50  element 
array  and  1  min  for  a  2 5 -element  array  are  anticipated 

It  is  possible  to  run  the  BASIC  program  in  its  uncompiled 
state  The  execution  of  the  program  w  ith  cache  memory  and 
the  fast  CPU  board  as  tne  only  enhancements  results  m 
computation  times  of  approximately  25  mm  for  a  50-element 
arras  and  4  5  mm  for  a  25  element  array 

A  copy  of  the  entire  program  is  available  from  the  authors 
Our  specific  implementation  in  HP  BASIC  utilizes  several 
hardware  and  software  devices  to  achieve  computational 
efficiency,  some  of  which  may  not  be  pertinent  to  other 
BASIC  operating  systems  running  on  comparable  machines 
Users  will  undoubtedly  find  it  necessary  to  make  modifications 
to  the  code  to  allow  it  to  run  on  other  HI’  equipment  or  in 
BASIC  tin  the  VAX 

V  Exampifs 

The  following  examples  demonstrate  the  utility  of  the 
current  algorithm  for  application  in  situ  and  provide  insight 
into  different  situations  that  might  arise  when  reshadmg 
cquispaced  arrays  with  failed  elements  If  optimal  reshading 
can  restore  the  array  's  original  design  sidelobc  level  by 
slightly  increasing  the  mainlobc  beamwtdth.  then  we  say  that 
the  optimal  reshadmg  has  been  effective  Optimal  re.shading  is 
effective  m  many  common  element  failure  situations  When 
the  arrav  is  severely  degraded,  optimal  reshadmg  is  less 
effective  but  is  stili  useful  in  reducing  the  negative  impact  of 
element  failures  These  examples  demonstrate  that  the  effec¬ 
tiveness  of  reshadmg  depends  upon  the  number  of  element 
failures,  as  well  as  the  location  of  the  failed  elements  within 
the  array 

Missing  elements  are  modeled  by  zeroing  the  appropriate 
weights  In  these  examples,  ,V  refers  to  the  number  of  intact 
array  elements.  M  is  the  number  of  beam  pattern  samples,  and 
A’o  is  calculated  by  using  the  equation  in  an  earlier  footnote 
Wc  define  the  mamlobe  width  to  be  twice  K0  rn  all  three 
examples 

A  Example  /.  Effective  Reshading 

This  example  demonstrates  that  reshadmg  can  restore  the 
original  sideiobe  level  of  an  array  response  by  slightly 
increasing  the  mamlobe  bcamwidth  In  a  25-element  equi 
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spaced  array,  originally  designed  tor  dB  Mdetobes. 

elements  2  and  4  have  failed  Iherelore.  V  23,  M  128. 
and  A„  -  0.68?**  Wc  first  keep  liu*  mainiobc  width  fixed  and 

allow  the  side  lobe  level  to  rise  See  Fig  1  !‘hc  peak  sidelobc 
level  has  risen  to  •  26.86  dB  below  the  mamlobe.  and  the 
main  lobe  width  is  unchanged  11  the  sideiobe  level  after 
reshadmg  is  too  high,  an  alternative  to  discarding  or  repairing 
the  array  is  to  broaden  the  mamlobe  bcamwidth  In  Fig  2.  A 
is  increased  to  0  775  and  the  peak  sideiobe  level  diminishes  t<» 

30  04  dB  below  the  mamlobe  A  trade  off  must  always  be 
made  between  an  enlarged  mamlobe  beam  width  and  an 
acceptable  peak  sideiobe  level  In  this  case  the  mamlobe  was 
increased  12  7  percent  in  order  to  recover  the  original  sideiobe 
level  Execution  times  on  the  HP9836C  are  between  1  and  2 
nun  for  Figs  1  and  2 

R.  Example  2‘  Moderated  Effective  Reshadmg 

This  example  is  taken  from  Siren  and  Nuttull  |1}  Because 
of  the  improvements  detailed  in  Section  HI.  above,  the  current 
algorithm  runs  faster  on  the  HP98  36C  than  on  the  VAX  !  I 
?80.  although  the  floating  point  multiply  time  on  the  HP  m  sis 
basic  configuration  is  roughly  200  times  slower  than  on  the 
V  AX 

Consider  a  linear  array  with  50  cquispaced  elements, 
•nmaily  designed  for  peak  -.idelobes  of  M)  dB  relative  to  the 
mamlobe  Fig  3  shows  the  classical  Dolph  Chebyshcv  beam 
pattern  with  30  dB  sidelobes  throughout  the  minimization 
range  |A(i.  (2*771  A  j  This  w.o  computed  uvng  the 
current  algorithm  m  6  1  !  miri  This  ideal  case  could  have 
been  computed  analytically  » 

Now  we  suppose  that  five  elements.  "'.22. 40.  a  U  SO.  of  the 
array  have  failed  The  optima:  response  after  reshadmg  the 
arras  is  shown  in  Fig  4  The  peak  sideiobe  level  has  risen  to 

25  51  dB.  but  we  have  maintained  mamlobe  bcamwidth  and 
retained  lull  steering  capability  In  this  example  A  4^  and 
M  128 

This  example  < Fig  4i  took  “4"  minutes  on  the  HP9K36C 
and  required  202  simplex  iterations  The  algorithm  of  Strcit 
and  Nuttali  required  38  4  mm  and  402  iterations  on  the  VAX 

Recovery  of  the  original  sidelobc  level  is  possible  t:ig  5) 
The  mamlobe  bcamwidth  must  be  increased  by  the  large  factor 
257  6  percent  (A..  -  0.87])  and  the  execution  of  this  task 
takes  8  98  min  and  requires  351  iterations  The  constraint  that 
all  the  weights  he  between  0  and  l  is  used  It  is  necessary  to 
Use  the  constraint  in  this  instance  because  otherwise  a 
dislocation  of  the  maximum  response  from  k  =  0  results  This 
dislocation  is  due  to  the  presence  of  loo  many  negatively 
weighted  elements 

C.  Example  3:  A  Severely  Degraded  Array 

This  example  shows  that,  for  severely  degraded  arrays, 
recovery  of  the  original  sideiobe  level  may  not  be  possible  by 
increasing  the  mamlobe  bcamwidth.  even  after  optimal  re 
shading  Consequently,  control  of  the  level  of  the  first  sideiobe 
must  be  relinquished  in  order  to  gain  control  of  the  level  of  the 
remaining  sidelobes 

Consider  a  25  element  array  with  elements  1 1  and  14  failed 
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ELEMENT  NUMBER 

Fig  I  Optimized  array  response  and  normalized  weights  for  25  elements 
with  elements  2  and  4  missing 
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ELEMENT  NUMBER 

Fig  2  Array  response  and  normalized  weights  for  Example  I  with  A'0  = 
0.775. 


The  original  sidelobe  level  is  -  30  dB.  Here  N  =  23,  M  =  minimization  interval  by  moving  A'o  far  enough  to  the  right 
128.  and  K0  =  0.6877.  Fig  6  shows  the  algorithm's  optimal  will  improve  the  approximation,  but  one  must  give  up  control 

response  to  this  configuration  11  is  a  significant  observation  of  the  first  sidelobe  to  reduce  the  others  to  acceptable  levels, 

that,  in  this  case,  small  perturbations  of  K0  will  not  affect  the  The  net  effect  of  losing  two  elements  so  close  to  the  center  is 
level  of  the  sidelobes  Only  when  the  first  sidelobe  is  that  negligible  emphasis  is  placed  on  the  remaining  center 

incorporated  into  the  mainlobe  beamwidth  ( K0  =  1 .27)  does  elements  (12  and  13)  and  the  rest  of  the  aperture  is  reshaded  as 

the  level  of  the  remaining  sidelobes  return  to  the  original  if  it  were  two  separate  arrays. 

desired  value  (see  Fig  7).  It  is  apparem  that  decreasing  the  This  situation  cannot  be  overcome  by  using  different 
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Fig  3  Classical  Dolph-Chebyshcv  array  response  and  normalized  weights 
for  ,V  -  50  and  -  30  dB  sidelobcs 
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Fig  4  Optimized  array  response  and  normalized  weights  for  50  elements 
with  elements  7.  22,  40,  43.  and  50  failed 


weights.  The  optimal  property  of  the  array  problem  formula¬ 
tion  and  solution  tells  us  that  no  weights  exist  which  can 
suppress  ail  the  sidelobes  below  a  certain  level  Thus  this  array 
has  lost  too  many  elements  and  performance  cannot  be 
restored  to  its  original  design  levels  merely  by  reshading. 

We  have  chosen  to  relinquish  control  of  the  first  sidelobe  to 


gain  cjntrol  of  the  level  of  the  remaining  sidelobes.  We  pick 
the  first  sidelobe  merely  for  ease  of  implementation;  modifica¬ 
tion  of  the  algorithm  to  forfeit  control  of  a  different  sidelobe 
could  also  have  been  done.  The  need  to  relinquish  control  of 
the  first  sidelobe  level  has  only  appeared  in  cases  of  severe 
array  degradation  due  to  element  losses. 
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Fig  5  Recovery  of  original  sidclobc  level.  Example  2  with  K0  -  0  8'1 


Fig  6  Optimal  array  response  and  normalized  weights  for  25  elements  with 
elements  li  and  14  failed 


VI  Conclusions 

Arrays  that  have  failed  elements  can  be  reshaded  to  obtain 
optimal  array  response  functions  Optimal  reshading  is  effec¬ 
tive  tn  many  common  element-failure  situations  When  the 
array  is  severely  degraded,  reshading  is  less  effective,  but  still 
can  be  used  to  reduce  the  negative  impact  of  element  failures 


Optimal  reshadtng  can  be  accomplished  in  situ,  quickly  and 
reliably,  on  portable  microcomputers  using  the  algorithm 
described  here.  Arrays  with  25  elements  routinely  run  in  less 
than  2  min  and  computation  time  for  a  50-element  array  is  less 
than  10  min  The  algorithm  can  be  applied  to  arrays  of  evenly 
or  unevenly  spaced  linear  geometry. 

The  above  examples  (and  others)  support  the  generally 
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Fig  7.  Recovery  of  original  sideiobe  level.  Example  3  with  K0  =  1.27. 


accepted  notion  that  failure  of  near-center  dements  is  more 
detrimental  to  the  airay  response  than  failure  of  near-edge 
elements. 

Another  application  of  Algorithm  635  is  to  arrays  of  planar 
and  arbitrary  three-dimensional  geometry.  Computation  times 
for  these  more  general  arrays  probably  will  depend  upon  N 
(number  of  sensors)  and  M  (  number  of  beam  pattern  samples) 
in  the  same  manner  as  for  linear  arrays. 
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Appendix  B 

DRIVER  PROGRAM  LISTING:  “Reshade” 
Lines  1  through  481 


1  OPTION  BRSE  1 

2  OUTPUT  2  USING  " tt ,  B "  ;  255 ,  75  1  CLEAR  SCREEN 

3  PRINTER  IS  CRT 

4  RAD 

c  REAL  Es  t  ore  k  50  > ,  IK  256  )  ,  Tbeami  256  1  .  Gr  r  e  =  (.  256  ’ ,  Hi  e n  ;  '  2 5C 

6  INTEGER  Io«xi t (10), It log(10>, Icsunt (5S), Ijtut  k  3f51 ' 

7  INTEGER  Ld  i  m ,  I ,  J ,  K ,  N ,  C5 ,  Somf  l  ag ,  Cac hf  1  ag ,  F  1  oat  f  1  ag,  H 

8  COM  /Prrs£/  Zradi i < 50) , Brad l t < 4 ) , Che b ( 1 6 ) , Z < 50 > , 2c ent r v SO  1 

9  COM  /firrssl/  Ref  <  256  )  ,  I  mf  <  256  )  ,  Reb  <  4 , 50  >  BUFFER , I  mb < 4 , 50 >  BUFFER  ,  Rebc  ent  r  < 
4 ) ,  I  mbc  ent  r  <  4  > 

10  COM  /Proj/  Bas  i  nwC  5 1 , 54 ) ,  Coss  i  n<  2 ,  1 025 ) ,  Rea<  256 , 50  )  BUFFER ,  I  mat  256 , 50  >  BUF 
FER,Cos45, Space 

11  COM  /'Par  am/  INTEGER  Nd > m, M , L , Logp , Nd imp  1 , Nd i mp4 , S 11 , Tme , M t s e I C 1 0 ■ 

12  CON  /Buffmult/  Co  1  real 50) , Col i ma< 50 ) , Co  1  re b < 50) , Co  1 t mb <50 ) 

13  COM  /Groups/  INTEGER  Nogroup, REAL  Sens  1 er , XgroupC 25 ) , Dgroup ,  Ik 50 > 

14  COM  /Groupsl/  Hydseni I 3 ] , Hydro$ C 3 3 

15  DIN  S 11  $  [  3  3  ,  Equ  t  *  [  3  3  ,  We  i  gh t  *  1 3  3 ,  Hegue t  *  1 3  3 ,  tie  uk  o*  C  3 3 ,  Bat  a_msus#C  18  3 ,  Ft  ler.a 
me$C103,Wgtstore*C33, Group_spac  e* C  3  3 

16  ! 

1?  Dat a_msus*= “ s I NTERNAL" 

18  Cac  hf 1 ag  =  0 

19  ON  ERROR  GOTO  24  !  POSSIBLE  ERRORS  IF  INTERFACE  NOT  PRESENT 

20  CONTROL  32 ,  1 ; 1  !  IF  CACHE  MEMORY  IS  PRESENT  IT  WILL  BE  UTILIZED 

21  OFF  ERROR 

22  STATUS  32, l;Stats 

23  IF  Stats  THEN  Cachflag=l 

24  Floatflag=0 

25  ON  ERROR  GOTO  Redo 

26  CONTROL  32,2|1  !  IF  FLOATING  POINT  CARD  PRESENT  IT  WILL  EE  UTILIZED 

2 7  OFF  ERROR 

28  STATUS  32 , 2 i St  at  s 

29  IF  Stats  THEN  Floatflag=l 

30  ! 

31  Redo  1  !  OBTAIN  INPUT  DATA 

32  LOOP 

33  IF  Ndt  m=0  THEN 

34  Nd i m=  1 6 

35  ELSE 

36  Nd i m=Nd t m+Tme 

37  END  IF 

38  REPEAT 

39  PRINT  "ENTER  TOTAL  NUMBER  OF  ELEMENTS/GROUPS  IN  ARRAY:  (3-50)  C  "t.VALt  ( 
Ndint)&“3" 

40  INPUTNdim 

41  UNTIL  Nd  t  m  >2  AND  Ndim<51 

42  i 

43  OUTPUT  2  USING  "  #  ,  B " ; 255 , 75  !  CLEAR  SCREEN 

44  REPEAT 

45  PRINT  "ENTER  NUMBER  OF  SENSORS  IN  EACH  GROUP:  C  "S,VALf<  Nogroup  >S, "  3  " 

46  INPUT  Nogroup 

47  UNTIL  Nogroup<26 

48  IF  Nogroup=0  THEN  Nogroup=l 

49  ! 

50  IF  Nogroup<>l  THEN 

51  REDIM  Xgroup(Nogroup) 

52  OUTPUT  2  USING  "#,B“;255,75  '  CLEAP  SCREEN 

53  REPEAT 

54  Group_spac e*= " " 

55  INPUT  “IS  ELEMENT  SPACING  WITHIN  THE  GROUP  CONSTANT7  <Y/N)  CY3",Grou 

p_spac  e* 

56  IF  LtN<Group_space*)=0  THEN 

57  Group_space$="Y" 

58  ELSE 

59  Group_space*=UPC*(Group_spacetC 1 3  > 

60  END  IF 

61  UNTIL  Group_space*=" Y”  OR  Group_space*="N" 
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63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85  ' 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 
101 
102 

103 
3" 

104 

105 

106  ! 

107 

108 

109 
(D'C) 

110 
1  1  1 
112 
1  13 

114 

115 

116 

117 

118  ! 

119 

120 
121 


IF  Gr oup_spac e*  =  " N “  THEN 
REPEAT 
H  =  0 

PRINT  "ENTER  POSITIONS  OF  SENSORS  IN  GROUP:" 

FOR  1=1  TO  Nogroup 

PRINT  "SENSOR  #  "  8.VAL*  (  I  >8.  "  :  " 

INPUT  Xgr oup  < 1 > 

IF  I>1  AND  Xgroupl I XXgroupC 1-1)  THEN  H=H*1 
NEXT  I 
UNTIL  H=0 
ELSE 
REPEAT 

PRINT  "ENTER  SPACING  BETWEEN  SENSORS  IN  CROUP:  C  "fcVALi  <.  I'gr  oup  .>  &  "  3  “ 

INPUT  Dgroup 
UNTIL  Dgroup>0 
FOR  1=1  TO  Nogroup 

Xgroup(I3=<I-l>*Dgroup 
NEXT  I 
END  IF 
ELSE 

MAT  Grr es=  <1.  ) 

END  IF 

OUTPUT  2  USING  "4, E"! 255, 75  !  CLEAR  SCREEN 

REPEAT 

Hydsen*= " " 

PRINT  "DO  YOU  WISH  TO  INCORPORATE  A  HYDROPHONE  SENSITIVITY'  <VrN>  C  N  3 " 
INPUT  Hydsen* 

IF  LEN<Hydsen*)=0  THEN 
Hydsen*="N“ 

ELSE 

Hydsen*=UPC*(Hydsen*C 1 3 ) 

END  IF 

UNTIL  Hydsen*="Y"  OR  Hydsen*="N" 

IF  Hydsen*="N"  THEN 
MAT  Hsens  =  (1,3 
ELSE 

OUTPUT  2  USING  "# , B " ; 255, 75  1  CLEAR  SCREEN 

REPEAT 

PRINT  "ENTER  THE  PHYSICAL  SENSOR  LENGTH:  <  METERS  )  C  "&-VAL  f  f  Sens  1  en  >  H  " 

I NPUT  Sens  1 en 
UNTIL  Sens  1 en >0 . 

REPEAT 

Hydro*= " " 

PRINT  "IS  HYDROPHONE  TO  EE  MODELED  AS  A  DIPOLE  OR  CONTINUOUS  SENSOR7 

CC3" 

INPUT  Hydro* 

IF  LEN ( Hydro*  3  =0  THEN 
Hydro*="C" 

ELSE 

Hydro*=UPC*<Hydro*C 1 3  ) 

END  IF 

UNTIL  Hydr o*= " C "  OR  Hydro*="D" 

END  IF 

OUTPUT  2  USING  "#,B";255,75  \  CLEAR  SCREEN 

REPEAT 

PRINT  "ENTER  TOTAL  NUMBER  OF  MISSING  ELEMENTSyGROUPS  I  "S.VAL*<  Tme  >8."  3 


122  INPUT  Tme 

123  UNTIL  Tme >  =  0  AND  Ndim-Tme>2 

124  ! 
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125  IF  Tme  >0  THEN 

126  REDIM  Miscl <Tme> 

12?  REPEAT 

128  OUTPUT  2  USING  "#,B";25S,?S  ■  CLEAR  SCREEN 

129  PRINT  “ENTER  MISSING  ELEMENT  GROUP  NUMBERS  (SEPARATED  EV  COMMA;' 

j  Mi se 1 <  *  > i "  3s” 

130  INPUT  Misel <*> 

131  MAT  SORT  Misel (*> 

132  H  =  0 

133  FOR  1=1  TO  Tme 

134  IF  Misel  <1X1  OR  Misti  CI)>Ndi«  THEN  H  =  H-t  i 

135  IF  I>1  THEN 

136  IF  Mi  sell  I  >=Misel < 1-1 )  THEN  H=H*1 

13?  END  IF 

138  NEXT  I 

139  UNTIL  H=0 

140  END  IF 

141  ! 

142  OUTPUT  2  USING  “ 4 . B " ; 255 , 75  I  CLEAR  SCREEN 

143  REPEAT 

144  INPUT  "ARE  ALL  ELEMENTS'  GROUPS  EQUISPACEB''  <Y/N)  IY]",Equi* 

145  IF  LEN ( Equ l * ) =0  THEN 

146  Equ i *= “ Y " 

14?  ELSE 

148  Equi*  =  UPC*<Equi*CU; 

149  END  IF 

150  UNTIL  Equ  t  *  =  " Y “  OR  £qui*=”N" 

151  ! 

152  REDIM  D(Ndim-Tme) 

153  New  ko:  ! 

154  ~ Symf 1 ag=6  1  FLAG  FOR  ARRAY  SYMMETRY 

155  IF  Equ l *= " Y "  THEN  '  EQUISPACED  ARRAY 

156  OUTPUT  2  USING  " « , B “ j 255 , 75  I  CLEAR  SCREEN 

157  IF  S1I =0  THEN  SI  1 =30 

158  REPEAT 

159  PRINT  “ENTER  ORIGINAL  SIDELOEE  LEVEL  <  D  B  >  i  <0  TO  501)  I - “ t VRL S C S 1 

“3" 

160  INPUT  SI  1  * 

161  IF  LEN  <  S 1  1  * ) <  >0  THEN  S 11 =ABS < VAL< S II  * ) ) 

162  UNTIL  SI  1  >-l  AND  SI  1 <5 1 

163  i 

164  OUTPUT  2  USING  "#,B"j 255,75  !  CLEAR  SCREEN 

165  REPEAT 

166  PRINT  “ENTER  ELEMENT  /-GROUP  SPACING  (METERS  >  <0-15>  t  "(kVAL*  (Space 
3  " 

167  INPUT  Space 

168  UNTIL  Space>0  AND  Space<=15 

169  ' 

170  N=Nd  irti'l 

171  R=10'(S!  1/-20)  '  CALCULATE  K0 

172  R2=R*R 

173  R3=SQR ( R2~  1 . > 

174  R5=(R+R3)A ( 1 . /N ) 

175  R6=(R-R3)/'(1./N) 

176  2o=(R5  +  R6>/'2. 

177  Ko=(2./'Space)*ACS(l./'2o) 

178  Kl=2.  *PI/-Space-Ko 

179  ! 

180  IF  Neuk o*= “ Y "  THEN 

181  OUTPUT  2  USING  "  It ,  B  "  ;  255 , 75  '  CLEAR  SCREEN 

182  REPEAT 

183  PRINT  "ENTER  K0:  C  "  &  VAL*<  Ko )  8,  ”  ]" 

184  PRINT  "SUGGESTED  VALUE  IS  : " ; PROUND ( Ko , -4 ) 

185  INPUT  Ko 

186  Kl=2. *PI/Space-Ko 

187  IF  Hydsen*=“Y"  OR  Nogroup>l  THEN 


I 


188 

189 

190 

191 

192 

193 

194 

195 

196 
19? 

198 

199 

200 
201 
202 

203 

204 
285 
206 
20? 
208 

209 

210 
211 
212 

213 

214 
?  1  * 
216 
21? 
218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 
24? 

248 

249 

250 

251 

252 

253 


PRINT  "ENTER  K  1 :  I  "  S.VRL*  <  k  1  “  1  “ 

INPUT  K 1 
END  IF 

UNTIL  Ko>0  AND  Ko<PI'Space  AND  Ko<Kl 
Nd i  m  =  Nd  i  m  +  Tnie 
END  IF 

C5  =  0 

FOR  1=1  TO  Ndim 
IF  T me >0  THEN 
FOR  J  =  1  TO  Tme 

IF  I=MiseKJ>  THEN  204  I 

NEXT  J 
END  IF 
C5=C5+1 

D<C5>=Space*<  I  -1  ) 

NEXT  I 

CALL  Symd < Nd  i  rn-Tme , Sy mf t ag , D < » ) > 

ELSE  . 

PRINT  "ENTER  ELEMENT 'GROUP  POSITIONS  (METERS  FROM  END /  j* 

PRINT  “SKIP  MISSING  ELEMENT'GROUP  POSITIONS." 

IF  Newko*=“Y“  THEN  Nd i  (n=Nd l  m  +  Tme 
FOR  1=1  TO  Ndim-Tme 
REPEAT 
H  =  0 

PRINT  "ELEMENT  'GROUP  "LVALTCDS."  C  "  S.VAL*  < D  <  I>H  ”  I  :  "  J 

INPUT  D<  I  >  I 

IF  I>1  THEN 

IF  t(I)<B<I-l)  THEN  H=H+ 1 
END  IF 
UNTIL  H=0 
PRINT  D ( I ) 

NEXT  I 

OUTPUT  2  USING  “#,B";255,?5  '  CLEAR  SCREEN 

REPEAT 

INPUT  “ENTER  KO:  ( RAD 'ME TER ) “ , Ko 
INPUT  “ENTER  K 1  *  (RAD'METER) " , K1 
UNTIL  Ki/no  mnD  Ko>=0 
CALL  SymdC Nd l m- Tme , Symf 1 ag , D ( * > 1 
END  IF 

Ndi m=Ndi m-Tme 
Nd i mpl =Nd i m  + l 
Nd i mp4=Ndi m  +  4 
IF  Eqv i t = “ Y “  THEN 
M  =  64 

C3=  <  K l -Ko ) '( 2. *M-1  . ) 

ELSE 
M  =  1 28 

C3=<K1-Ko)'(M-1. > 

END  IF 
Logp=5 

PI  =  <2'-Logp>  +  l 

OUTPUT  2  USING  "#,B“j255,75  (  CLEAR  SCREEN 

REPEAT 

Negwet$=“  " 

INPUT  “WILL  YOU  ALLOW  NEGATIVE  WEIGHTS  ?  <Y'N>  C Y 1 " , Negwet t 
IF  LEN<Negwet*)=0  THEN 
Negwet  *= " Y “ 

ELSE 

Negwet»=UPC*<Negwet*[ll) 

END  IF 

UNTIL  Neguet*=“Y“  OR  Negwet*="N" 

IF  Neguet*  =  "Y"  THEN  !  NO  CONSTRAINTS:  <WjX  =  l 
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25? 
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260 
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312 
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314 
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ELSE  1  1  CONSTRAINT:  <  SUM(  W  j  >  -  .  5 1  (  = .  5 ,  <  W  .1  -  .  5  )  =.5 

L=  1  1  CHANGE  L  AND  REDIM  APPROPRIATE  ARRAYS 

END  IF  !  FOR  MORE  CONSTRAINTS 

Ld  i  m  =  Mf;X  (  1  ,L) 

'  REDIMENSION  INPUT  ARRAYS 

i 

REDIM  I  oex  i  t  <  Logp  1 ,  I  j  suit  ( 3,  Ndi  mp  1 1 ,  1 1  I  og  <  Logp  1 ,  I  c  cunt  <  Nd  i  (n  ’ ,  Zr  ad  i  i  <  Nd  t  (« 
REDIM  Brad  i  i  < Ldi  m  1 ,  CMeb (Logp  1 ,  Est  ore < Nd i  m  1 ,  Z (Nd i  m  1 ,  Zcent  r  <  Nd  i  fn  1 
REDIM  Basinv(Ndimpl,N.-inp41,Cossin(2,Pl>>Rea(M,Ndim>fIfina(M,Ndirti> 

REDIM  Ref  (Ml , I rnf CM) , Reb(Ldi m, Ndi ml , Imb(Ldim,Ndif«l,F:ebcentr<Ldi(n> 

REDIM  I  mbc  entr(Ldim),  D<  Nd  i  m  1  ,  U(  M  1  ,  Tbeam  (  2*M  >  ,  Grres  (Ml ,  Hsens  (Ml 
REDIM  Colrea(Ndiml,Col  i  ma(  Nd  i  ml  ,  Co  1  reb(  Nd  i  iri ) ,  Co  I  i mb(Hd i ml 

i 

MAT  Bas i nu  =  (0.1  •  INITIALIZE  COMMONS 

MAT  Coss in=  (0.1 
MAT  Z  =  <0.  ) 

MAT  Cheb=  (0.  1 
MAT  Co  1  r ea=  (0.1 
MAT  Colima=  (0.) 

MAT  Co  1 r eb  =  (0.) 

MAT  Col  i  rnb  =  <0.  1 

I 

IF  Neguiet$  =  "Y"  THEN 
MAT  Reb=  (0. > 

MAT  I  mb  =  <0.) 

MAT  Rebcentr=  (0.1 
MAT  Imbcentr=  <0.1 
MAT  Brad  1 i =  (0.1 
MAT  Zcentr=  <0. ) 

MAT  Zradi 1=  (1. > 

ELSE 

MAT  fieb  =  (1. ) 

MAT  I mb=  (0.) 

MAT  Rebcentr=  (.5) 

MAT  Imbcentr=  (0.) 

MAT  Bradi 1 =  < .  51 
MAT  Zcentr=  (.51 
MAT  Zradi i =  ( .  5) 

END  IF 

FOR  1=1  TO  M 

U(I)=Ko+C3*(I-ll 
NEXT  I 

| 

IF  Hydsen*=“Y"  THEN 
MAT  Hsens=  <0. ) 

IF  Hydro*= " D "  THEN 
FOR  J=1  TO  M 

Const  = . 5+ . 5*C0S(U<  Jl*Sens 1 en 1 
Hsens(Jl=Const»Const 
Const=-.5*SIN(U(Jl*Senslenl 
Hsens(Jl=SQR(Hsens(Jl+Const*Constl 
NEXT  J 

ELSE  !  CONTINUOUS  SENSITIVITY 

FOR  J=1  TO  M 

Hsens(Jl=ABS(3IN(U(Jl*SensIen/2. >/(U(JHSenslen/2.  11 
NEXT  J 
END  IF 
END  IF 

f 

IF  NogroupOl  THEN  !  CALCULATE  GROUP  RESPONSE 

MAT  Grres=  <0. > 

FOR  J=1  TO  M 
Gr i m=0 . 

FOR  1=1  TO  Nogroup 


!  GENERATE  U  ARRAY 

i  CALCULATE  SENSITIVITY  TERM 
!  DIPOLE  SENSITIVITY 
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385 


Grres<  J  >  =  Gt-resK  J  >  +  COS  (  U<  J  .<  *Xgroup  <  1  >  > 

Gr  i  ra  =  Gr  i  m-SIN<U<  J  >*Xgroup<  I  >  > 

NEXT  1 

Grres  (  J  >=SQR  (  Gr  res  <  J  >  *Gr  res  (  J  )  +  Gr  !  fn*Gr  i  m  >  'Nr.  gr  oup 

NEXT  J 
END  IF 

FOR  J  = 1  TO  M 

Ref  <  J  )  =C0S  <  D(  Nd  l  m  )  >  i  GENERATE  F  ARRAY 

Imf<J)=-SIN<D<Ndtm)*lKJ)> 

FOR  1=1  TO  Ndim 

ReavJ, n=ReftJ)-COSiD(l)*UCJ)>  '  GENERATE  HR  ARRAY 
Iioa<  J,  n  =  imf  ( J>+SINCD<  I  >«U(  J)  i 
Rea<  J , I ) =Rea( J , I )*Grres<J)*Hsens(J> 

I  mat  J ,  I)=Ima<J,  I>*Grres<J)*Hseris(J> 

NEXT  I 
NEXT  J 

FOR  1=1  TO  M 

Ref<Il=Ref<Il«Gi'res<I>*H=eris<I> 

Imf<I)  =  Imf(I)*Grres<I)*Hseris<I> 

NEXT  I 


N=Nd l m- 1 

Itlog<l)=20*N  !  MAX  ITERATION  COUNT 

I oex i t  < 1 ) =0  IPRINT  OPTION 

Ts=T I MEDATE  (INITIALIZE  TIME 

CALL  Kapro»<N,Itlog<*),IoexltC*),l  jsvil  (•*  1 1 
Te=TIMEDATE-Ts  (EXECUTION  TIME 

Estore<N>=Cheb( Logp ) 

I  count  <N )  =  I 1 1 og< Logp ) 

Zsum=0‘  !  CALCULATE  FINAL  UEIGHT=1-SUM  OF  ALL  OTHERS 

Zsum  =  SUM  <  Z ) 

Z  <  Nd i m )  =  1 .  -Zsum 

IF  Symflag  THEN  i  SYMMETRIZE  WEIGHTS 

FOR  1  =  1  TO  I  NT  < ( Nd i m ) r2  > 

P=(Z(I )  +  Z<Ndi m-I  +  1 )  >r2 
Z<  I)=P 

Z ( Nd i m- I + 1 ) =p 
NEXT  I 
END  IF 

IF  Tme  >0  THEN 

REDIM  Z C Nd i m+ Tme ) 

FOR  1=1  TO  Tie 

FOR  J  =  Nd  i  m+ 1  TO  HiseK!)  STEP  -1 
IF  J  >M i se 1 <I>  THEN  Z<J)=Z<J-1) 

NEXT  J 

ZCMi  seU  I )  )=0. 

NEXT  I 
END  IF 

IF  Equ i t  =  "  Y  "  THEN 

CALL  Cal cbeam<Ndim,M, Tme, Mi sel <»), Space, Z < * >  ,  TbeanK*  >  > 

ELSE 

CALL  UnsymbeamfNdi m, M, Tme, Mi se l<»>,Ko,Kl,Z<*>,Tbeam<*>> 

END  IF 

Zinax  =  MAX<  Z  <  *  1  ) 

IF  Zmax<>0  THEN 

FOR  1=1  TO  Ndim+Tme  !  NORMALIZE  WEIGHTS  TO  1 

ZCI >=Z(I)xZmax 
NEXT  I 
END  IF 


* 


fti 


» 


» 


» 


i 


B-8 


3S6  1 

387  CALL  We  1  ght _pt  1  ot  <, Nd  1  m  +  Tme ,  M-'2 , 8 ,  Tins  ,  M  i  se  1  '  *  1 ,  Spac  e  ,  Ko ,  (  1 ,  C  3 , 3  ■  -  .  T  te  at.-  -  * 

) , Equ l *  > 

3S8  OUTPUT  2  USING  " # , B" ; 255 , 75 

389  PRINTER  IS  PRT 

390  PRINT  USING  "0" 

391  BUMP  GRAPHICS 

392  ! 

393  CONTROL  CRT, 12j2 

394  FOR  1=0  TO  9 

395  ON  KEY  I  LABEL  “ “  GOSUB  Dummy 

396  NEXT  I 

397  ON  KEY  1  LABEL  '  CONTINUE  "  GOTO  Comp 

398  ON  KEY  2  LABEL  "  KEYS  OFF-'ON  "  GOSUB  Flip_key 

399  LOOP 

400  END  LOOP 

401  Dummy:  ! 

402  RETURN 

403  F 1 i P_k  ey :  ! 

404  Kef  1 ip=<Kef1  ip+1 >  MOD  2 

405  IF  Kef  1  ip  THEN 

406  CONTROL  CRT, 12; 1 

407  ELSE 

408  CONTROL  CRT, 12; 2 

409  END  IF 

410  RETURN 

4 1 1  Comp :  ! 

412  GRAPHICS  OFF 

413  OFF  KEY 

414  We i ght *  =  "  " 

415  REPEAT 

416  INPUT  “WOULD  YOU  LIKE  A  LIST  OF  THE  WEIGHTS''  <V-N,  C  Y I "  ,  We  i  ght  * 

417  IF  LEN  <  We i ght  t  >  =0  THEN 

418  We i ght  *  =  " Y " 

419  ELSE 

428  We i ght *  =  UPC*< We i ght *11 3 > 

421  END  IF 

422  UNTIL  We l ght  *= " Y "  OR  We i ght *  = " N" 

423  ! 

4  24  CALL  Pr  i  ht  i  hput  s  ( Cac  hf  1  ag ,  F  1  oat  f  1  ag ,  Nd  i  m ,  Tme  ,  S 1  1  ,  1 1  1  og '  Logp ,  M  i  s  e  1  '»  1 , 

Spac  e , Ko , K 1 , Te , Cheb  t  *  > , 2  <  *  3 , Equ i * , We i ght  * , Ne  gwe  t  *  > 

425  ' 

426  OUTPUT  2  USING  "#,B";255,75  >  CLEAR  SCREEN 

427  Neuk  o*  =  "  11 

428  REPEAT 

429  INPUT  "WOULD  YOU  LIKE  TO  CALCULATE  A  NEW  ko  or  K1  TO  GIVE  A  DIFF.  DEAN 
WIDTH  < YyN>  CN3",Neuko* 

430  IF  LEN  <  Neuk  o*  >  =  0  THEN 

431  Neuk  o*  =  " N " 

432  ELSE 

433  Neuko*=UPC*<Neuko*C 1 3 ) 

434  END  IF 

435  UNTIL  Neuko*="Y”  OR  Neuko*="N" 

436  IF  Neuk  o*= " Y"  THEN  GOTO  Neu_ko 

437  i 

438  Wgt s t or e *  =  "  “ 

439  REPEAT 

448  INPUT  "WOULD  YOU  LIKE  TO  STORE  THE  WEIGHTS  IN  A  DATA  FILE''  (Y/N>  I N  3  "  , 

Wgt  st  ore* 

441  IF  LEN(Wgtstore*3=0  THEN 

442  Wgt s t ore*= "N " 

443  ELSE 

444  Wgt s t or e*=UPC* CUgt st ore* t  1 3 > 

445  END  IF 

446  UNTIL  Wgt  s t  ore*= " Y “  OR  Wg t s t or e*  =  “N" 

447  IF  Wgtstore*=”Y"  THEN 
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H 


448 

444 

450 

451 

452 

453 

454 

455 

456 
45? 

458 

459 

460 

461 

462 

463 

464 

465 

466 
46? 

468 

469 

470 

471 
4  72 

473 

474 

475 

476 
47? 
47s 

479 

480 

481 


OUTPUT  2  USING  ' # , B" ; 255 , 75  1  UEf.R  SCREEN 

INPUT  “ENTER  FILENAME  FOP  WEIGHT  FILE:  i  1 0  C  HuFm1'  TEE- 
OUTPUT  2  USING  "  ** ,  B";255,  75  1  CLEAR  SCREEN 

INPUT  "ENTER  MASS  STORAGE  DEVICE:  I  :  INTERNAL  J  '  ,  I  =,»  i 
IF  Equi*=”Y"  THEN 
IF  Nd : m  >38  THEN 

CREATE  BDAT  F i 1 enamet&Dat a  msust, 2, 256 
ELSE 

CREATE  BDAT  Fi  1  e  n  am  e  t  S,  D  a  t  a  msust. 1,256 
END  IF 

ASSIGN  0Stordat  TO  F  i  1  eriarnc  tS.Dat  a_msust 
OUTPUT  @St  ordat  jNdim  +  Tfne,Spac£,2».*> 

ELSE 

SELECT  Nd i m 
CASE  >4? 

CREATE  BDAT  Fi  ienimellDsta  msust, 4,256 
CASE  >31 

CREATE  BDAT  F  i  1  ename*8,Dat  a  msust, 3, 256 
CASE  >15 

CREATE  BDAT  Fi lenamet&Data  msust, 2, 256 
CASE  ELSE 

CREATE  BDAT  F i 1 enamet&Dat a  msust, 1,256 
END  SELECT 

ASSIGN  GStordat  TO  F  i  1  er.ame  t«,Bat  a_r»s ust 
OUTPUT  8Stordat;Nd\m+Tme,D(*>,Zc*> 

END  IF 

ASSIGN  9Stordat  TO  * 


END  IF 
GRAPHICS  ON 
PAUSE 

GRAPHICS  OFF 
END  LOOP 
END 


RETURNS  TO  Redo  AT  PROGRAM  BEGINNING 


» 
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